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^ ■ Abstract 



Influence of uniaxial small-scale anisotropy and compressibility on the stability of scaling regime 
and on the anomalous scaling of structure functions of a scalar field is investigated in the model 
Vi*/ I of a passive scalar field advected by the compressible Gaussian strongly anisotropic velocity field 

■ with the covariance oc 5{t — i')|x — x'p^ by using the field theoretic renormalization group and 
the operator product expansion. The inertial-range stability of the corresponding scaling regime 
is established. The anomalous scaling of the single-time structure functions is studied and the 
corresponding anomalous exponents are calculated. Their dependence on the compressibility pa- 
rameter and anisotropy parameters is analyzed. It is shown that the presence of compressibility 
leads to the decreasing of the critical dimensions of the important composite operators, i.e., the 
anomalous scaling is more pronounced in the compressible systems. This result is demonstrated 

', for the structure function of the third order. All calculations are done to the first order in e. 

\0 ■ Introduction 
O 

■ The so-called "rapid change model" of a scalar field passively advected by a self-similar Gaussian 
5— correlated in time velocity field introduced by Kraichnan [Ij and number of its extensions play 
the central role in the theoretical investigation of intermittency and anomalous scaling [21, fS]. The 

■ main reason for this is the experimental fact that the deviations from the statements of the famous 
^ . classical Kolmogorov-Obukhov phenomenological theory (see, e.g., [21 [3]) are more noticeable for 

^ I simpler models of passively advected scalar quantity (scalar field) than for the velocity field itself (see 
[4] and references cited therein) and, at the same time, the problem of the passive advection of a scalar 
field is much easier from theoretical point of view than the original problem of anomalous scaling in 
the framework of the Navier-Stokes velocity field. Within the rapid change model of a passive scalar 
advection the systematic analysis of the corresponding anomalous exponents was done for the first 
time on the microscopic level. For example, in the so-called "zero- mode approach" to the rapid change 
model (see survey paper ^ and references cited therein) the anomalous exponents are found from the 
homogenous solutions (zero modes) of the closed equations for the single-time correlations. 

One of the most effective approach for studying self-similar scaling behavior is the method of the 
field theoretic renormalization group (RG) [U [6] . This method can be also used in the theory of fully 
developed turbulence and related problems, e.g., in the problem of a passive scalar advection by the 
turbulent environment (see [El [71 E] and references cited therein) . 

In [9] the field theoretic RG and operator-product expansion (OPE) were used in the systematic 
investigation of the rapid-change model of a passive scalar. It was shown that within the field theoretic 
approach the anomalous scaling is related to the very existence of the so-called " dangerous" composite 
operators with negative critical dimensions in the OPE (see, e.g., Refs. [6l|8] for details). 
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In subsequent papers a few generalizations of the rapid-change model towards more realistic ones 
were done. For example, in [10] the field theoretic RG and OPE were applied to the rapid change 
model of a passive scalar advected by the (5-correlated, Gaussian strongly anisotropic velocity field. It 
was shown that in the small-scale anisotropic case the anomalous exponents of the structure functions 
and correlation functions are nonuniversal and they are functions of the parameters of anisotropy, 
Besides, they form the hierarchy with the leading exponent related to the most "isotropic" operator. 
On the other hand, in [TT] the influence of compressibility and large-scale anisotropy on the anomalous 
scaling behavior was studied in the aforementioned model and the anomalous exponents of higher-order 
correlation functions were calculated as functions of the parameter of compressibility. It was shown 
that the presence of small-scale anisotropy is more pronounced for larger values of the compressibility 
parameter a. Prom this point of view, compressible systems are very interesting to be studied. 

In the end, in [l2j the combined effects of the small-scale anisotropy and compressibility on the 
anomalous scaling of the structure functions of a passive scalar within the rapid-change model were 
studied by using the simplest possible generalization of the incompressible and anisotropic tensor 
structure of the velocity field correlator to the compressible one where compressibility was introduced 
only to the isotropic part of the anisotropic tensor structure. In present paper, we shall investigate 
more general case of the model where compressibility will be introduced through all anisotropic tensor 
structures of the velocity field statistics. Pirst of all, we shall establish stability of the scaling regime 
of the model and coordinates of the corresponding infrared (IR) stable fixed point will be found as 
functions of the compressibility and anisotropy parameters. These results will be then used in the 
analysis of the asymptotic behavior of the single-time structure functions of a passively advected scalar 
field. 



Formulation of the model 

We shall investigate the model of the advection of a passive "tracer" Q{x) = 9{t,x.) which is described 
by the following stochastic equation 

dt9 = uoAO - {vidi)e + /, (1) 

where dt = d/dt, di = d/dxi, A = 0^ is the Laplace operator, vq is the coefficient of molecular 
diffusivity (a subscript will denote bare parameters of unrenormalized theory), Vi = Vi{x) is the i-th 
component of the compressible velocity field v(x), and / = f{x) is a Gaussian random noise with zero 
mean and correlation function 

Df = {f{x)f{x')) = 5it - t')C{v/L), r = X - x', (2) 

where parentheses (...) hereafter denote average over corresponding statistical ensemble. The concrete 
form of the noise defined in ([2|) will not be essential in what follows. The only condition which must 
be satisfied by the function C(r/L) is that it must decrease rapidly for r = \r\^ L, where L denotes 
an integral scale related to the stirring. 

In real problems the velocity field v(x) satisfies Navier-Stokes equation but, in what follows, we 
shall work with a simplified model where we suppose that the velocity field obeys a Gaussian statistics 
with zero mean and pair correlation function 

{vAx)v,{x'))^Dl{x-x') = D,5{t-t') j (^(^^^^%^-p[zk(x-x')], (3) 

where d is the dimension of the space, k is the wave vector, and Dq is an amplitude factor related to 
the coupling constant of the model (expansion parameter in the perturbation theory, see, e.g., [12] ) 
by the relation Dq/vq = = A^^, where A is the characteristic UV momentum scale. The parameter 
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of the energy spectrum of the velocity field < e < 1 is taken in such a way that its "Kolmogorov" 
value (the value which corresponds to the Kolmogorov scaling of the velocity correlation function in 
developed turbulence) is e = 2/3, and 1/m is another integral scale. In our uniaxial anisotropic and 
compressible case the tensor Rij(k) is taken in the following way \10\ I12j 



<y.2nsni{Pis + aQis){Pji + aQji) , (4) 

where we denote ^/j = n • k//c, i-ij(k) = 6ij — kikj/k'^ is common isotropic transverse projector, 
Qij = kikj/k'^ is the longitudinal projector, the unit vector n determines the distinguished direction 
of uniaxial anisotropy, a > is a free parameter of compressibility, and ai, are parameters char- 
acterizing the anisotropy. From the positiveness of the correlation tensor one finds restrictions on 
the values of the above parameters, namely, a\^i > — 1- 

The stochastic problem ([TJ-Q can be rewritten in a field theoretic form with the following action 
functional [SJ [6] 

= J dt d'^x 6' [-dt - Vidi + I/O A + xoi^oln • df'\ ^-\Jdti d'^xi dt2 d'^xa 

(i;i(tl,Xi)[L>j;.(fi,Xi;t2,X2)]"^t;i(t2,X2) -0'(tl,Xi)L>^(tl,Xi;t2,X2)6''(t2,X2)) , (5) 

where D^j and d{j are given in ([3]) and ([2]) respectively, 9' is an auxiliary scalar field (see, e.g., [6]), 
and the required summations over the vector indices are implied. In action ^ the term with new 
parameter xo is related to the presence of small-scale anisotropy and the introduction of this term 
is necessary to make the model multiplicatively renormalizable. Model ([5]) corresponds to a stan- 
dard Feynman diagrammatic technique (see, e.g., [l2] for details) and the standard analysis of 
canonical dimensions shows which one-irreducible Green functions can possess UV superficial diver- 
gences. The functional formulation ([5]) gives possibility to extract large-scale asymptotic behavior of 
the correlation functions after an appropriate renormalization procedure which is needed to remove 
the UV-divergences. 

Influence of anisotropy and compressibility on the scaling regime of 
the model and on the anomalous scaling 

The IR scaling regimes of the model are given by the IR stable fixed points of the corresponding RG 
equations [6l[7l[8]. The fixed points of the RG equations can be determined from the requirement that 
all the so-called beta functions of the model vanish and the IR stability of the fixed point is given by 
the requirement that all the eigenvalues of the matrix of the first derivatives Qij = dPi/dCk, where 
Pi denotes the full set of beta functions and is the full set of charges of the model, must have 
positive real parts. In our case the coordinates of the fixed points are given by the following system 
of equations 

l3g{g*,X*,aj,a,d,£) = g^{-2£ + -fl) = 0, (3^{g^,x*,aj,a,d,£) = x*(7i -72) = 0) (6) 

for J = 1, 2 and variables with star denote the fixed point values of the corresponding variables. The 
explicit form of the so-called gamma functions 7^, i = 1,2 is as follows 

71 = 1^1 , 72 = ^^2, (7) 
2 2x 

where we denote g = gS(i/{2TT)'^ , = 27r'^/^/r(d/2) is the d— dimensional sphere, and Ai^2 have the 
following explicit form 

. (d + 2)(d - 1) + ai(d+ 1) + 02 + a(d + 2 + 01-202) + 0^02 

= d(dT2) ' 
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Figure 1: The restriction on the values of the anisotropy parameters ai and 02 for different values of 
the compressibility parameter a for space dimensions d = 2 and d = 3 given by the physical condition 
X* > —1- The figures show that the allowed regions of the parameters ai and 02 decrease when the 
compressibility parameter a increases. 



-2(qi + 02) + d'^a2 + 2a{ai + ^02) + 2a^a2 

= d{dT2) • 

Thus, possible fixed points are found as solutions of the system of algebraic equations ^ and their 
IR stability is determined by the positive real parts of the eigenvalues of the matrix 17 = {^lik}- Using 
the explicit expressions for gamma functions d?]) we have found the coordinates of the nontrivial fixed 
point (trivial fixed point with g^, = is not interesting for us) 

Hd + 2)e 

^* {d + 2){d-l)+ai{d+l) + a2 + a{d + 2 + ai-2a2) + a^a2' 
_ -2{ai + 02) + d'^a2 + 2a{ai + da2) + 2a'^a2 

~ {d + 2){d - 1) + ai{d + 1) + 02 + a{d + 2 + ai - 2a2) + a^a2 ' 

At fixed point the matrix of the first derivatives has the eigenvalues Ai,2 = 2e. It means that the 
point is IR stable if e > 0. It is also the only condition to have > (of course, together with the 
physical assumptions: ai_2 > — 1 a a > 0). Besides, the fixed point value of x must satisfy physical 
condition x* > ~1- This fact leads to restrictions on parametrical space of the model, namely, on 
parameters of anisotropy ai_2 and compressibility parameter a. For space dimensions d = 2 and d = 3 
these restrictions are shown in Fig.[TJ The existence of these restrictions is not present in the model 
with simpler definition of the anisotropy and compressibility which was studied in [12j. 

Further we shall use the above obtained results to determine the anomalous scaling of the single- 
time structure functions of the scalar field inside the so-called inertial interval which is defined by 
inequalities [3l[6l [lOl [12] . They are defined as follows 

S^(r)^([^(t,x)-0(t,x')]^), (12) 

where denotes the order of the structure function and r = |r| = |x — x'|. We shall not discuss details 
of the corresponding analysis. It will be given elsewhere but we only stress that the anomalous scaling 
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Figure 2: Dependence of the critical dimension A[3, V\/e on the anisotropy parameters ai and 02 for 
different values of the compressibility parameter a and for d = 2>. 

is given by the existence in the model of the so-called dangerous composite operator with negative 
critical dimensions within the OPE [3 [3 [8]. In our case the main contribution to the anomalous 
scaling is given by the operators F[N,p\ = ■ ■ ■ di^O {diOdiOY with N = 2n + p (see, e.g., [TOl [12]). 

After rather long additional renormalization procedure [101 US] one comes to the following final 
expression for the inertial range behavior of the single-time structure functions 

S^{v) = Z)-^/V^(i-) T.{Cn'.. (r/L)^[^''^l + ...}, (13) 

N'<N P 

where A[A^,p] denote the critical dimensions of the operators F[N,p] (more precisely, they are eigen- 
values of the corresponding matrix of the critical dimensions, see, e.g., [E]), p obtains all possible 
values for given N' , Cn'<p are numerical coefficients which are functions of the parameters of the model, 
and dots means contributions by the operators others than F[N,p] (see, e.g., [HI [TO] for details). 

Our aim is to analyze the influence of compressibility on critical dimensions of anisotropic operators 
and to find the answer on the question whether the compressibility makes the anomalous behavior 
more pronounced or not, i.e., whether the presence of compressibility decreases the critical dimensions 
of corresponding operators or not. 

In Fig. [21 the dependence of the critical dimension of the principal eigenvalue (related to the 
isotropic case, see p!0]) of the corresponding matrix of critical dimensions of the composite operators 
-F[3,p] is shown which play the central role in the asymptotic behavior of the structure function 5*3 (it 
is well known that the structure function 5*2 has not anomalous behavior |10l I12j . therefore we start 
the analysis from 5*3). It can be seen that in the case of the structure function 5*3 the anomalous 
scaling is present (the existence of negative values for A [3, 1]) and the effect is more pronounced for 
more compressible system. Detail analysis of the model for the structure functions of higher order will 
be given elsewhere. 
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Conclusions 

In present paper we have studied the influence of uniaxial small-scale anisotropy and compressibility 
on the anomalous scaling of the structure functions of a scalar field in the framework of the Kraichnan 
model by using the field theoretic RG and OPE. Using the RG technique we have shown the existence 
of the scaling regime within the inertial range which is defined by the corresponding IR stable fixed 
point of the RG equations. We have found restrictions on the anisotropy parameters which are 
related to the compressibility of the system. The restrictions become larger when compressibility 
parameter increases. Besides, using the OPE we have found the asymptotic form of the structure 
functions and, using the concrete example of the third order structure function, we have shown that 
the compressibility of the system makes the effects of anomalous scaling more pronounced, i.e., the 
critical dimensions of the corresponding composite operators are smaller in compressible case than in 
the incompressible case. Thus, we can conclude that the compressible environment is more suitable 
for experimental study of the anomalous scaling of the structure functions of a passive scalar than 
incompressible systems. 
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